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In this paper we study the critical behavior of an A'^-component c/j'^-model in hyperbolic space, 
which serves as a model of uniform frustration. We find that this model exhibits a second-order phase 
transition with an unusual magnetization texture that results from the lack of global parallelism 
in hyperbolic space. Angular defects occur on length scales comparable to the radius of curvature. 

This phase transition is governed by a new strong curvature fixed point that obeys scaling below 
the upper critical dimension d^c = 4. The exponents of this fixed point are given by the leading 
order terms of the 1/N expansion. In distinction to fiat space no order 1/A^ corrections occur. We 
conclude that the description of many-particle systems in hyperbolic space is a promising avenue to 
investigate uniform frustration and non-trivial critical behavior within one theoretical approach. 


I. INTRODUCTION 


Field theory and statistical mechanics in geometries 
with negative curvature are of increasing interest. While 
a direct application to the spacetime of our universe 
seems to require a positive cosmological constant, a 
wide range of many-particle problems are closely tied to 
problems with negative spatial curvature. For example, 
field theories in hyperbolic space are increasingly studied 
because of its direct relation to anti-de Sitter space. 
The latter is essential for the duality between strong 
coupling limits of certain quantum field theories and 
higher-dimensional gravity theories[IH3]. The scaling 
behavior near critical points in hyperbolic space, being 
the Wick-rotated version of anti-de Sitter space, may 
therefore be of relevance in the analysis of strong 
coupling theories. On the other hand, networks like 
the Bethe lattice, that have been studied early on in 
the statistical mechanics of phase transitions HD and 
that have received renewed interest in the context of the 
dynamical mean-field theory of correlated fermions [71 [S], 
quantum spin glasses jS], or bosons [TU], can be considered 
as a regular tiling of the hyperbolic plane HU. To be 
precise, if one considers a regular tiling {p, q}, where p 
refers to the degree of a polygon and q to the number 
of such polygons around each vertex, then the Bethe 
lattice with coordination number q corresponds to 
{cx),( 7 }. All regular tilings of the hyperbolic plane with 
{p — 2){q — 2) > 4 are possible [TT]. Obviously the 
square lattice {4,4}, the triangular lattice {3,6}, and 
the honeycomb lattice {6, 3}, i.e. the only possible tilings 
with regular polygons of the two-dimensional fiat space, 
are just excluded. It was already stressed in Ref. m 
that hyperbolic tiling might be used to interpolate 
between the mean-field behavior of the Bethe lattice and 
a lattice that might be close to the square or honeycomb 
lattice. This may offer an alternative approach to study 
corrections beyond dynamical mean-field theory. A tiling 
of three-dimensional hyperbolic space with dodecahedra 
is shown in Figure (see |12|1. 

Finally, effects of uniform frustration are often captured 


in terms of certain background gauge fields or by embed¬ 
ding a theory in curved space m [M] . An interesting 
case of tunable uniform frustration is found by studying 
a given fiat-space problem in curved space with inverse 
radius of curvature k, an idea that was introduced in 
|15fil8| . Here, the problem of packing identical discs 



Figure 1: Tessellation of three-dimensional hyperbolic space 
by dodecahedra (created with |12| 1 

was studied in a hyperbolic plane. While in fiat space 
packing in hexagonal close-packed order is possible, in 
hyperbolic space this order is frustrated by the fact 
that gaps open up between neighboring discs. This 
facilitated the study of packing properties as a function 
of frustration. The latter can be varied by changing the 
spatial curvature k. Thus one might capture packings 
that are not allowed in fiat space, as it occurs for clusters 
with icosahedral local order, in terms of a non-frustrated 
model that is embedded in a curved geometry. These 
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ideas were also employed in studies of glass transitions in 
hyperbolic space (imiiD]), where the authors performed 
molecular dynamics simulations on the hyperbolic plane 
for a Lennard-Jones liquid and found that the fragility 
of the resulting glass is tunable by varying k. 

Given these applications of negatively-curved geometries, 
it is an interesting question to ask how phase transitions 
of classical and quantum models will behave in such 
curved spaces. Significant numerical work has been 
devoted to studies of classical spin models in hyperbolic 
space. The thermodynamic properties of Ising spins 
placed on the vertices of lattices in hyperbolic space were 
studied in |21fEl] . In order to perform Monte Carlo 
simulations on finite two-dimensional lattices a nega¬ 
tively curved background is created by tessellating the 
hyperbolic plane with regular n-gons. All these works 
have found the phase transition to follow mean-field 
behaviour. In particular, various critical exponents were 
measured and found to numerically coincide closely with 
mean-field exponents. One should, however, keep in 
mind that the detailed protocol for measuring the critical 
exponents in these works is somewhat different from 
the usual flat space protocol. The mean-field behavior 
is supported by a Ginzburg criterion for 0^-theories 
in hyperbolic space that was discussed in Ref. |25| . 
There this behavior was rationalized by arguing that the 
Hausdorff dimension of hyperbolic space is infinite. 

The problem now arises to address the question of phase 
transitions in three-dimensional hyperbolic space. In 
particular, it is natural to ask, whether there is scaling 
as in flat d-dimensional space, or if the transition is 
genuinely of mean-field type. If phase transitions in hy¬ 
perbolic space were of mean-field nature, it would imply 
for systems that are below their upper critical dimension 
in flat space, that an arbitrarily small curvature k would 
lead to a violation of scaling. An alternative possibility 
is that scaling continues to be valid in hyperbolic space 
with a new fixed point characterized by new exponents. 
The above numerics would in that case indicate that 
some exponents take their mean-field values. Because 
of hyperscaling this cannot be the case for all critical 
exponents. In the case of a new fixed point there are 
obvious questions: what is the universality class and 
what are the critical exponents? 

We answer these questions in this paper by study¬ 
ing analytically the problem of an A^-component 
continuum ^^-theory in three-dimensional hyperbolic 
space in a large-A expansion. In the discussion section 
we comment on the generalization to different dimen¬ 
sions. Section of this paper contains the exposition 
of the 4>‘^ model in hyperbolic space. We find that the 
theory possesses a second order phase transition, that 
scaling is obeyed below the upper critical dimension 
of the flat space and that the exponents are given by 
the leading order terms of the 1/A expansion. To be 
specific, we find that the leading order 1/A-corrections 
to the exponents vanish. In addition we give general 


arguments that support the conjecture that all higher 
order 1/A-corrections should vanish as well. The 
technical steps of our calculation are as follows: Using 
the momentum space analysis of section EYi we locate 
the critical temperature of this phase transition and 
find the magnetization texture of the ordered phase. 
In Section |III| we discuss the character of the phase 
transition and present the results of the calculations of 
the exponents rj, v and 7 to lowest order in 1/A. In 
order to meaningfully identify the critical point of this 
model, it is convenient to formulate the problems in 
momentum space. As this representation in hyperbolic 
space does not seem to exist in the condensed matter 
literature, we develop the necessary parts in section IV. 
With this formalism it is now possible to deal with the 
order 1/A correction to the critical exponents 77 and 7. 
We find that the exponents rj^v and 7 at lowest order 
are those of three-dimensional flat space and not those 
of a mean-field transition. However, in distinction to 
flat space, the 1/A corrections vanish. The absence of 
higher-order corrections is found to be the consequence 
of the finite curvature of hyperbolic space, which ex¬ 
ponentially cuts off fluctuations of wavelengths longer 
than the curvature radius. This is in agreement with the 
general remarks on the regulating behavior of hyperbolic 
space by Gallan and Wilczek in |26| . 

The critical exponents satisfy scaling and we discuss in 
the final section how our results can be understood from 
the scaling of the free energy in the presence of finite 
spatial curvature. 

As a further result we calculated the magnetization 
texture of the ordered state of this model. We find that 
uniform magnetization develops in regions of size 1/k. 
Due to the lack of the concept of global parallelism in 
hyperbolic space (CHIIITI), these regions will necessarily 
be uncorrelated in their magnetization direction. 

II. MODEL AND BACKGROUND GEOMETRY 

The model we are considering is an A-component 
theory given by the action 

+ j ( 1 ) 

with summation over i, fj, and v implied. Thus we are con¬ 
sidering here the three-dimensional version of (/)^-theory. 
Generalizations to different dimensions are straightfor¬ 
ward, as we discuss in the final section. In the action, 
is the metric of three-dimensional hyperbolic space, 
which is a maximally symmetric space with negative cur¬ 
vature, characterized by a single parameter, the curva¬ 
ture K. The quantity g is the metric determinant and 
assures the proper transformation property of the action. 
Hyperbolic space can be defined as one of the two (equiv¬ 
alent) simply-connected three-dimensional manifolds of 


3 


points satisfying 

2 , 2,2 2 1 

Xi -\- X2 + — -^ ( 2 ) 

Kj 

inside Minkowski space. The coordinates Xi are the carte¬ 
sian coordinates of Minkowski space. To derive a more 
convenient formulation, the points may be parametrized 

by 


Xi = 

— sinh KT sin 9 cos (p 

K 

(3) 

X2 = 

— sinh KT sin 9 sin p 

(4) 

2^3 = 

— sinh KT cos 9 

K 

(5) 

X4 = 

— cosh KT. 

(6) 


Minkowski space has a metric that is given by 

ds^ = dx\ + dx 2 + dx\ — dx\. (7) 


This formula follows from the fact that the vectors are 
parallel-transported such that the angle between the 
geodesic curve and the vector is a constant. Since hyper¬ 
bolic triangles have angles which sum to tt — ('|28|h 

we are left with the angular defect stated in ([^. Thus 
the ordered state in hyperbolic space will in general be 
more complicated. Inside the radius of curvature, where 
A < 1/k^, a uniform direction may be meaningfully de¬ 
fined. 

In order to determine the nature of the ordered state, we 
study the symmetry-broken state of the action S at the 
lowest order in a 1/A^ expansion, i.e. we peform a saddle 
point analysis of the partition function 

Z = [ D4)exp{—S), (10) 


and then include higher-order fluctuations in a system¬ 
atic fashion. Here, the action S is given by: 


S = /3 


^^21 


cf) {^IQ - cf) + ^ (cf) ■ cf)f] . ( 11 ) 


This induces an intrinsic metric on the hyperbolic space 
with line-element 


= dr^ + ^ sinh^ kt {dO^ + sin^ . (8) 

In the limit k —>■ 0 we regain three-dimensional flat space. 
An additional length scale, 1/k, is present in hyperbolic 
space, which is ultimately responsible for the non-trivial 
magnetization texture that we derive below. Note that 
our results can straightforwardly be applied to quantum 
phase transitions in hyperbolic space, if one of the spa¬ 
tial coordinates is considered as imaginary time after the 
usual Wick rotation. 


III. PHASE TRANSITION AND 
MAGNETIZATION TEXTURE 

In three-dimensional flat space the model that we con¬ 
sider is known to possess a second order phase transi¬ 
tion, where the ordered state corresponds to a symmetry- 
broken phase with uniaxial magnetization. Geometri¬ 
cally this is not possible in hyperbolic space (IIHIET]), 
since here a global direction is not a well-defined concept. 
Consider, as shown in Fig. three locally magnetized 
patches A, H, C, which are the corners of a hyperbolic 
planar triangle. The direction of the order-parameter 
at A may be parallel-transported to B and C along the 
geodesics AB and AC, respectively. If now we continue 
the parallel-transport from B to C along BC, the two 
transported magnetization directions will not match. In¬ 
stead, there will be an angular defect 9 between the two 
magnetizations that is proportional to the enclosed area 
A of the hyperbolic triangle: 


0 = Ak^ 


(9) 


We rewrite this by performing a Hubbard-Stratonovich 
decoupling of the {4> ■ <p)^ term, whereupon the action 
becomes 



dV 


cj) (/io -k iX{x) 




Now we integrate out all (j)i with the exception of the one 
component, along which the spins near a chosen point 
order and which we will label a{x). Moreover, we in¬ 
troduce a source held h{x) for a{x). This leads to the 



Figure 2: Three magnetized patches in a plane, forming a 
hyperbolic triangle. It is not possible to give a meaningful 
definition of the direction of magnetization, since parallel- 
transport of a vector from A to B and A to C, will not result 
in vectors that will match upon parallel-transport from B to 
C or vice versa. 
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action 


S = 13 dV 


a{x) (/^o + *A(a;) - V^) a{x) + ^A^(£c) 


+ [\og{fio + iX{x) - V^)] 


—/? J dV h{x)a{x) 


for the partition function Z = J Z^crDA exp(—S'). 

The saddle point solutions are determined by the condi¬ 
tions = 0 and . = 0 , which result in the two 

d(T{x) dA(a5) ’ 

equations 

{n{x) — V^) cr(a;) = h{x) (12) 

+ ^x’(x) ( 13 ) 

with fj.(x) = fj-o + iXix). 

We will use these equations to work out the critical ex¬ 
ponents including higher-order corrections in section IV. 
Here, we will only analyze equation (IT^ to find the sus¬ 
ceptibility x{x,x') = 5a{x)/5h{x'%^^,)^Q, which, by 
virtue of ( [I^ , satisfies 

(/r(£c) - V^) x{x, x') = 5{x - x'). (14) 

As we approach the phase transition from high temper¬ 
atures, /i(a:) may be assumed to be homogeneous. Us¬ 
ing the formalism of the next section, this equation may 
be transformed into momentum space whereupon it be¬ 
comes 


Xp = 


1 

/r -f _|_ j,2 ’ 


(15) 


where p > 0 are the eigenvalues of the Laplacian. Note 
the presence of the ‘mass term’ in the denominator. 
This is a consequence of the fact that the Laplace oper¬ 
ator in hyperbolic space has a gapped eigenvalue spec¬ 
trum, as will be seen below. A criterion for the presence 
of the phase transition is the condition that Xp diverge. 
The highest value of p when this happens is p = —k^, 
where the p = 0 mode of the susceptibility diverges. Thus 
the phase transition takes place at p = —k^ with an order 
that is determined by the p — 0 Fourier mode. In con¬ 
trast to fiat space, the p = 0 eigenmode of the Laplace 
operator cannot be one of homogeneous order, in agree¬ 
ment with the foregoing argument about angular defects. 
Instead, it corresponds to a diminishing of the magneti¬ 
zation a along the one direction, that we chose not to 
integrate out. In other words, due to the lack of a global 
direction of magnetization, focussing on one component 
of the Wcomponent vector, entails that one is eventu¬ 
ally considering projections of the magnetization vector 
instead of the full vector. The diminishing of this projec¬ 
tion takes place according to the formula 


cr(r) = (To 


nr 

sinh nr ’ 


(16) 


where r is the geodesic distance from the origin, where 
the unintegrated component and local magnetization di¬ 
rection coincide and (Tq is the magnitude of the magne¬ 
tization at the origin. The phase transition corresponds 
to the formation of infinitely many patches, more pre¬ 
cisely three-dimensional regions, of characteristic sizes 
1 /k, which have nearly uniform magnetization. The de¬ 
cay of (j{r) in Eq. (16 1 does not imply a decay of the 
magnitude of the order parameter, but must be inter¬ 
preted as the order parameter rotating away from the 
chosen direction of the vector (p. 


IV. MOMENTUM SPACE REPRESENTATION 


We come now to the technical part of this paper that 


will allow us to analyze the saddle point equations ( 121 , 


(131 and compute critical exponents. In order to make 


progress with the calculations, it is convenient to obtain 
the momentum space representation of functions that are 
translationally invariant in hyperbolic space. Let '^{dpQ) 
be a given function of the geodesic distance between two 
points P and Q. The functional dependence on the two 
points will not have an arbitrary form, but will rather 
be expressed through the geodesic distance dpq between 
these two points. This distance is the length of the 
geodesic curve connecting these points. Explicit com¬ 
putation of this length yields the formula 

cosh ndpQ = cosh kt cosh nr' — sinh nr sinh nr' cos 7 
cos 7 = COS0COS0'-I-sin0sin0'cos ((/) — 0'). (17) 


The fact that such a function ip depends on the six co¬ 
ordinates not in an arbitrary way, but only through the 
geodesic distance, allows us to expand ipidpq) in terms 
of the eigenstates of the Laplace operator in hyperbolic 
space. Since hyperbolic space may be defined as the set 
of all points equidistant from the origin in Minkowski 
space, this Laplace operator is identical to the one ob¬ 
tained by writing down the 4-dimensional Laplace oper¬ 
ator in angular coordinates and restricting the distance 
from the origin to a constant. A similar situation was 
considered by Fock [29] , who studied the problem on a 3- 
sphere embedded in 4-dimensional euclidean space. The 
eigenfunctions of the Laplacian on this 3-sphere are the 
generalized spherical harmonics of three angles. Their 
full description was given in [29]. We find the eigenfunc¬ 
tions of the Laplacian in hyperbolic space by multiplying 
one of the angles in Fock’s solution by the imaginary unit, 
a prescription sketched briefly in an appendix of m- 
The hyperbolic Laplacian is given by 

1 

A = -^— drisinh^ Krdrip)-\ -^-As 2 

sinh^Kr ^ ^ sinh^iKr) 

As 2 = -^de (sindelp) + ^\-dlip. (18) 

smtl sin 0 ^ 

The eigenfunctions are then given by 

Ipqlm {r, 0,(p) = {Kr)Yim {0, (p) (19) 
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with eigenvalues 


radial part of the eigenvalue equation 


Atpqim{r, 0, (j)) = -( k ^ + q‘^)tpqlm{r, 9, 4>)- ( 20 ) 

Here the 1)^ ^.re the ordinary spherical harmonics on the 
2 -sphere and the Hp/ are special functions that solve the 


„ , d ^ 1(1 + 1)k^„ , > ,9 9,„ , 

—^Hgi -I-2KC0th«;r—Hqi-- Ilqi{K.r) = -(k^ + q^)Ilqi{K.r). 

d+ dr sinh"^ Kr 


( 21 ) 


r 


The solutions can be expressed in a Rayleigh-type for¬ 
mula 




sinh'a; / 


Ml \d(cosha;)^+^ 


where 


(-)' 

(-Y + 1^ 


(-Y+f 


\K/ 


\K/ 


cos{px) (22) 


(23) 


is a normalization constant. The differential equation 
being of the Sturm-Liouville form, these functions satisfy 
the orthogonality relation 

OO 

J dr sinh^ (kt) n,/(Kr)n,/i(Kr) = ^S{q - q'). (24) 


A. Addition theorem 


The eigenstates 'ippim(r,d,(j>) satisfy an addition the¬ 
orem, which was derived for the 3-sphere by Fock |29j . 
In the latter case of the 3-sphere this formula is fully 
analogous to the addition theorem for two-dimensional 
spherical harmonics. Again, by multiplying one of the 
angles by the imaginary unit, we obtain the correspond¬ 
ing addition theorem valid in hyperbolic space 


= ^£( 2 i + i)n„Mn„(K/)p,(cos7), 


(25) 


where the Pi{cosj) are the Legendre polynomials in 
cosy. 

As a demonstration of the use of this formula, let us de¬ 
rive the magnetization texture of the p = 0 eigenmode 
given in (16l. The eigenbasis expansion of x reads 


X(r,r')= f dp^i2l+ l)xpIlpi(Kr)Ilpi(Kr')Pi(cosj). 
i 

(26) 

Insertion of Xp from ( |15[ ) into this equation at the critical 
point p, = yields the real-space form of y. Now we 
construct the real-space form of only the p = 0 mode, 
which gives 

a{d{r,r')) = \iuiy^ (21 + l)^Upi(Kr)Ilpi(Kr')Pi(cos^) 

p—>0 p^ 


= CTq 


Kd(r, r') 
sinh Kd{r, r') ’ 


(27) 


as claimed. 


B. Extraction of coefficients and inversion formula 


The identity (251 will be crucial in obtaining the expan¬ 


sion coefficients of a given function ip(d() of the geodesic 
distance. This distance being a non-negative quantity, 
the value of ij) for negative arguments is irrelevant. In 
particular we may redefine tp for negative arguments such 
that it becomes an even function. This allows us to 
Fourier expand ip as follows 


+ 00 


'0(d) sinh Kd = 


27r^ 


dp 0pPsin(pd), (28) 


where we have chosen to split off a factor of np/ (27r^) in 
the definition of the expansion coefficient for later conve- 
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nience. Inserting (25) we obtain 


tjj{d) 


+00 


2tt^ sinh nd 


dp ippP siii{pd) 



— OO 


y^(2Z + l)'ippIlpi{Kr)Ilpi{Kr')Pi{cos'y) 
1=0 


(29) 


and have thereby managed to expand the arbitrary func¬ 
tion Ip in the new basis with coefficients 


+ 00 

V'p = — /" dxipdx]) smh{Kx)e~^^^. (30) 

up J 


Let us briefly comment on the structure of the expan¬ 


sion. Note that in (29) the expansion coefficient ipp has 


no dependence on In fact, the statement of (29) is 


that any function that depends on the set of coordinates 
{r,9,(p), {r',9',(p') only through the geodesic distance of 
the two points, can have no explicit I or m dependence 
of ipp. 

Conversely, to find the real-space function ip{d) from the 
knowledge of the coefficients ipp in the expansion (29) we 
use (28) which results in the inversion formula 


V'(d) 


ti 

27r^z sinh hid 


+00 

J dp^ppppp^. 


— OO 


(31) 


C. Convolution theorem 


Let / and g be two-point functions that depend on the 
geodesic distances dpg and dga, respectively. When we 
multiply these functions and integrate Q over all of hy¬ 
perbolic space, the resulting function h can only depend 
on the geodesic distance between points P and R. This 
convolution will in general be difficult to carry out in real- 
space. The fact that the 11^;^ and the Yim are orthogonal 
functions, however, allows us to reduce the convolution 
of / and g to a multiplication in momentum-space. This 
is seen explicitly by rewriting the relation 


dVgf{dpg)g{dgii) = h{dpp) 


(32) 


in the momentum representation (29) with expansion 


coefficients fp,gp,hp and using the orthogonality rela¬ 
tions for the radial functions and the spherical harmonics. 
Then this convolution formula translates into 


fp9p — ^p- 


(33) 


The solution of the Dyson equation below will require 


r 


knowledge about the momentum-space representation of 
the Dirac 5-function in hyperbolic space, which we denote 
by either Spg or S{r,r'). We define this function by 
the condition that convolution of an arbitrary function 
ip(dpg) with 6{dgp) must yield tp{dpp). Translating this 
condition into momentum space, we immediately read off 
from (33) the relation 6p = 1 and obtain thereby 


+ 00 


^PQ = ^ / dpJ2i‘^l + ^) 


1=0 


X Ilpi{Kr)Ilpi {nr') Pi (cos 7 ). 


(34) 


Conversely, however, the multiplication of two functions 
in real-space does not translate into a simple convolu¬ 
tion integral in momentum-space, but rather a double¬ 
integral. Given the product 

f{dpg)g{dpg) = h{dpg) (35) 


the corresponding momentum-space equation is found by 
employing the representation (31) 


OO OO 


hk = J J ‘^dfp9qPq^a.nh 


p + q — k ^ 
2k 


(36) 


— OO —OO 


i.e. instead of a single integral a double integral with 
kernel is obtained. 

In the limit fc —>■ 0 the symmetry properties of fp and g^ 
may be used to rewrite this kernel as 


lim hk 
k=,0 


OO 



— OO 



P9fp9q 

cosh^ ■ 


(37) 


This formula will be used in section |VI| in evaluating the 
corrections to the critical exponent 7 . 


V. CRITICAL EXPONENTS 


The previous formalism may now be employed to an¬ 
alyze the saddle point equations (12) and (13). These 


equations describe the physics of the model at lowest or¬ 
der in 1/N. Given this fact, all exponents that we derive 
below represent the lowest order contribution in 1/N. In 
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principle these are only the lowest order terms in an ex¬ 
pansion of the critical exponents in a power series in 1/N. 
In flat three-dimensional space there are indeed further 
corrections. The main result of this paper, derived in 
section lYll is to establish the absence of such corrections 
for r] and 7 in three-dimensional hyperbolic space. 

We begin by computing the bare Green’s function (p) 
in momentum space for constant p in the action for a(x). 
The real-space definition of the bare is obtained 

by inverting the action, i.e. 

{p-V^)G^°\r,r') =S{r,r'). (38) 

Inserting 


B. Critical Exponents v and 7 


We now study the behavior of the correlation-length 
as the critical temperature is approached from the dis¬ 
ordered regime by examining the saddle point equation 
(13l. In the regime p > —k^, the magnetization will be 


zero. We may therefore set cr = 0 and obtain 


M = Mo + 


— g' 

2Tr^ J + + 

0 


uT uT ! - 


(43) 


G(°)(r,r') = / dpY^{2l + \) 

X G^°'> ip)Ilpi {iir)Ilpi (nr')Pi (cos 7 ) (39) 
and the representation of ^(r, r') in (34), it is found that 

= (40) 


{fj, + +p2) 


At the critical point, p = — we have a power-law 
dependence on p. Employing the inversion formula, we 
And the real-space dependence 


G(°)(r,r') 


K r pPP'^ 

27r2isinhK(i J ^ {fj, + + p"^) 

— 00 


^ li+K?d 

27r sinh nd ’ 


(41) 


where d is the geodesic distance between r and r'. Evi¬ 
dently, even at the critical point the Green’s function de¬ 
cays exponentially, in accord with the previous remarks 
about parallel-transport. 


A. Critical Exponent 7 

Let us now proceed to the evaluation of the exponent 
77. At the critical point pL = the power-law form 
of the curved-space bare Green’s function in ( [40| agrees 
with the flat-space limit. The exponent 7 may therefore 
be defined through the relation 

G(P) « (42) 

We see that the bare G^°^(p), i.e. the lowest order form 
of the Green’s function in an 1/A^ expansion, has rj = 0. 


At T = Tc, where p = —k^, we have 
2 . 


(44) 


which allows us to remove po from (43). Now close to 
p > —K^, the quantity y/p + k? dominates over p + 
Thus we neglect the latter and And 




I T-T, 

TT Tn 


A. 


(45) 


The length ^ that diverges at the critical point is defined 
by 


e = 




(46) 


which is the natural length scale of the problem as follows 
from the correlation function Eq.([40| in its eigenbasis. At 
the same ti me, the real-space form of the correlation func¬ 
tion in Eq.(41l shows that G^°^(d) decays exponentially 
beyond the curvature 1 /k, even for = 0 . However, 
what matters is the eigenbasis of the correlation function, 
where G becomes scale invariant at Tc- The decay of the 
correlation function in real-space will, however, have im¬ 
pact on the 1/7V corrections of the critical exponents. 
Defining a critical exponent v through we And from 
(451 and (|4^ the exponent 


e-(T-T,)-i^72=l. 


(47) 


According to ( [15| the zero-momentum susceptibility is 

1 


X<j=o = 


p + 


(48) 


and by using (451 we And 7 = 2 . 


We emphasize here explicitly the fact that both expo¬ 
nents are not mean-field exponents. The latter are given 
by vmf = 5 and 7 mf = 1 - 



V 

7 

K = 0 

V + 0 


2-|l^+0(^) 

K 7^ 0 

0 


2 


Table I: Critical exponents in flat and hyperbolic space. 


VI. CORRECTIONS TO CRITICAL 
EXPONENTS 

Corrections to the critical exponents 7 and 7 are found 
by inspecting the self-energy. In flat space this calcula¬ 
tion is described in m and we And that the general 



























procedure carries over to hyperbolic space. This proce¬ 
dure consists in first determining how a correction to a 
critical exponent would manifest itself in the self-energy 
and calculating the order 0{1/N) diagrams to see if such 


contributions are present. 

We start with the action (111. 
real-space reads 


The Dyson equation in 


G{rp,rR) = G^°'>{rp,rR) + J drq j 


and is converted to 


G-i(p) = G(°)-i(p) + E(p) (50) 

by application of the convolution theorem. We follow m 
in rewriting this equation as 

G-i(p) = + I](0)) + E(p) - E(0) (51) 

and redefining the new inverse bare Green’s function 
G°“^(p) to be the first term, i.e. 

G^ ^ (p) = /io -|- E(0) = ji . (52) 

This has the advantage that at T = Tc the ‘mass’ fi + 
of the bare propagator vanishes. As a consequence of this 
redefinition, self-energy insertions in diagrams now take 
the form E(p) — E( 0 ) instead of E(p). 

The large-structure of the model allows us to restrict 
ourselves to a small number of diagrams. The calculated 
corrections will be exact to order 1/iV. The coupling con¬ 
stant u in the action is multiplied by a factor 1/A^. Due 
to the presence of N fields, there is a summation over 
the field index at every {(j> ■ ^)^ interaction vertex. We 
represent this interaction term by a dashed line. On the 
other hand, a summation over the field index at every 
vertex produces a factor N. Thus the series of bubbles 
connected by —u/N interaction lines, as shown in Figure 
are all of order 1/7V and need to be included for con¬ 
sistency. We denote this sum by a wiggly line and use 
the symbol D{p). It satisfies the relation 


D{p) 


u/N _ 2 

i + fn(p) '^"ivn(p) 


(53) 


where n(p) is the polarization operator. The approxi¬ 
mation comes from the large-u limit, which we will be 
considering from here on. In real-space 11 is given by 


2-v/ fi+K^d 

n(r,r') = G^ryf = (54) 

471"^ smh nd 

where d is the geodesic distance between r and r'. Using 
(301 we find 




^ I 2 2k k 


(55) 


drq, G^°\rp,rQ)Y,{rQ,rQ>)G{rQ,,rR) (49) 


where i^{z) = ^ logr(z) is the digamma function. With 
(531 we find 


Dp{p) = 


2ttp 


N 


Im 






(56) 


We present the calculation of the correction to rj in de¬ 
tail. The calculation of the correction to 7 is much more 
tedious and is only sketched. The result in both cases 
is the absence of any corrections due to the regularizing 
character of finite curvature. 


A. Order 0{1/N) correction of 7 


As we have seen 77 = 0 at lowest order in 1 /N. We now 
determine the 1/A^ correction to this result. We have 
defined 77 in ( [4^ . Such a correction would manifest itself 
in the self-energy. For large-iV this critical exponent can 
be expanded and reads 


G(p) = 


A-’' 


A-^ 


At: p"^ »? 47 r — 77 logp 


(57) 


i.e. a correction would lead to a p^ logp term in the self¬ 
energy and could be found as the coefficient of such a 
term. In fiat space there is indeed such a correction. We 
will now show that this p^ logp term of fiat space is reg- 
ular ized in hyperbolic space. The polarization operator 
(551 at the critical point p = —k? becomes 


n(p) = lto„h(^p). 


(58) 


In fiat-space the p^ log p contribution is produced by the 
diagram in Figure |^. To write down this term we take 
the real-space form of D{p), which is obtained from (531 
with (30) 


D{dpQ) = 


8k'^ cosh{KdpQ) 
TVtt^ sinh^(KdpQ) 


(59) 


wvw 


+ ---■ 



Figure 3: Dyson equation for the screened interaction D{p) 
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The self-energy in real-space is obtained by multiplying 
D with the bare Green’s function 


E(dpQ) = D{dpQ)G^°\dpQ) (60) 

_ 8 k ^ cosh(pdpQ) 

Ntt^ sinh®(KdpQ) 


Using again formula 
turn space as 


(30), we can write this in momen- 


S(p) 


8 (k^ +p'^) 
3iV7r2 



(62) 

y/ypy 51 
+p2 ^ ' 2j ’ 


where 7 is the Euler-Mascheroni constant. Taking the 
flat-space limit At —>■ 0 with fixed p, we obtain the asymp¬ 
totic relation 




The appearance of k in this formula is owed to the fact 
that we measure all momenta in units of At. 

In the opposite regime, where p tends to 0 for fixed cur¬ 
vature, we have instead of a logarithmic divergence the 
finite value 


This regularizing behavior of the finite curvature is 
shown in Figure where we defined a quantity r]{p) = 

, which in flat space would yield a finite 


d 


d log p 


S(p)-S( 0 ) 
p- 


r]. In hyperbolic space at sufficiently small p, i.e. long 
length-scales, the log behavior of the self-energy (62) is 
cut off and p{p) is suppressed to 0 . 


the full Green’s function 

G-'(m(T)) = po{T) + K^ + J:ipiT)) 

= Ai(r)-kAtU (65) 

In the following all calculations will be for zero external 
momenta, thus we have suppressed the momentum argu¬ 
ments. Subtracting from this equation the same equation 
evaluated at T = Eg, we And 

piT) + At" - E (^(T)) + E (0) = poiT - T,). ( 66 ) 

We have found 7 = 2 at lowest order in section |V B[ 
Writing I /7 = 1/2 — A and using the definition of the 

exponent via po{T — Tc) ~ {p{T) + for T near T^, 

we obtain 

- E (AA(r)) + E (0) ^ + At2 (67) 

-Aa/ p{T) + At2 log {p{T) + At") 

valid near T >Tc. The first term on the right-hand side 
is an 0(1) term and was already obtained in section [V| 
It is produced by a diagram, which is obtained from di¬ 
agram & by removing the internal wiggly line. 

There are two 0(1/A) diagrams that have to be consid¬ 
ered in computing the correction A, shown in Figures^ 
and|^. In flat space it can be shown that the diagram in 
Figure only gives a ~ ^ log p correction, whereas the 
diagram in Figure]^ in fact yields a finite A. We shall 
see now that in hyperbolic space neither diagram yields 
a contribution to A, as both logarithmic divergences are 
regularized by k. 


1. Diagram (a) 


We begin with the diagram in Figure [^. We denote 
this self-energy part by Ea(/r). Using eq. (37) we And 


B. Order 0{1/N) correction of 7 

The exponent 7 is found from the divergence of the 
susceptibility at p = 0. We And at zero-momentum for 



Figure 4: Relevant 0{1/N) diagrams. In flat space, diagram 
(a) contributes to 7 and diagram (b) to 7 . 


Nri(p) 

0.25 
0.20 
0.15 
0.10 
0.05 
0.00 

Figure 5: Here we have deflned a function ri{p) = 

j' s(pKS( 0 ) j^ ^ whereas in hy¬ 

perbolic space it is always regularized by curvature and tends 
to 0 at long length-scales. 
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A A 


Ea(M) - Sa(0) = ^ j dq j dp[Gp{^^)Dq{^^)-Gp(0)Dg{0)]—^^-J^——. 

A A ' 2k. d ) 


-A -A 


( 68 ) 


The hyperbolic cosine effectively cuts off contributions 
with \p — q\ ^ K. Hence the g-integral may be approxi¬ 
mated by limiting the integration to this region. Insertion 
of Gp{p) and Dq{p) from (401 and (56 1 and subsequent 
expansion in p results in a large number of elementary 
integrals. All logarithmic terms stemming from these in¬ 
tegrals are of the form log(cK^ -I- p), where c is either 1 
or 4. In other words, for finite k no terms proportional 
to \ogp are present. 


2. Diagram (b) 


Similarly the y/p log p divergence of the diagram 
in flat-space is regularized. The self-energy expression 
corresponding to this diagram is obtained by noticing 
that diagram & is the result of attaching to the self¬ 
energy in & two legs of the interaction vertex. It is 
correspondingly given by 


Sh(/r) = ^ qq' Gqip)"^ [T,aiq\p) 


1 


cosh' 


2 TTiq-q') 






[{q' — 2 kI'kY -\- p + -I- 2kI'kY -\- p + k^] 


Pa(g',M) - Sa(0,Ai)] 


(69) 


where the same approximation as before has been made. 
In the integrand the momentum-dependent self-energy 


Ea(g,/i) is required. According to (36 1 this is given by 


A A 

Sa(g', p) = J dp' J dp [Gp{p)D'p{p) - Gp(0)Dp(0)] pp' tanh ( j ^ 

-A -A 


Inside the p-integral we approximate the tanh-function in 
the region \p + p' — q\ < — hy its argument and outside 
this region by the sign-function. Then the p-integral may 


be carried out without a cutoff and we are left with a p'- 
integral. Insertion of (701 into (691 and integration over 
q' results in 


^b{p) 



^ (g'+ 2 a)2 + 4A2 


-(- R{k) . 


(71) 


where R{k) denotes terms that tend to 0 with a —>■ 0. 
The first term on the right-hand side reproduces for k = 
0 fully the flat space formula for E;,. This self-energy 
contains th e ^f p log p term. For finite k, however, the 
integral in (711 is fully regularized and a ^//2 log p term is 


avoided. Thus, we conclude that no singular correction 
to the p dependence of the self-energy emerges, i.e. the 
exponent 7 is also unchanged compared to the leading 
order 1/N expression given above. 
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VII. DISCUSSION 

The aim of this paper was to investigate critical phe¬ 
nomena in hyperbolic space. Our key finding is that for 
a ^^-model embedded in hyperbolic space a new fixed 
point emerges at finite curvature k. If k > 0 the critical 
exponents are governed by the strong curvature limit. 
Interestingly, these exponents are given by leading or¬ 
der terms of the 1/iV expansion. Thus, while the nu¬ 
merical values of the exponents are now simpler, they 
continue to obey hyperscaling below the upper critical 
dimension. The physical state in the symmetry-broken 
regime is characterized by an unusual magnetization tex¬ 
ture. This texture consists of regions of size of the order 
of the radius of curvature 1/k where the vector (p has 
nearly uniform direction. Beyond this region the finite 
value of the curvature starts to play an important role, 
since a global direction in hyperbolic space is not a well- 
defined concept. It is therefore not possible to establish 
a uniform direction of the magnetization vector. In fact, 
as we have illustrated in Figure the parallel transport 
of a local direction from region AtoB and then from B 
to C is not the same as the direct parallel transport from 
A to C. It is this lack of transitivity which is the origin 
of the resulting magnetization texture. 

The fact that the k ^ 0 values of exponents are differ¬ 
ent from the fiat space k = 0 limit may be understood 
using standard crossover arguments as we now show us¬ 
ing the example of magnetic susceptibility. Let /(f, /i, n) 
be the singular part of the free energy density, where 
t (X {T — Tc)/Tc measures the distance to the critical 
point and h is the external field. Then the following scal¬ 
ing transformation holds 

/(t, h, k) = b~^f b^^h, bn^ , (72) 

with exponents Vf for the correlation length and scaling 
dimension of the conjugate field yf = PfSf/vf that refer 
to the fiat space (k = 0) limit. The curvature is a relevant 
perturbation with positive scaling dimension, i.e. the in¬ 
frared behavior is governed by the infinite curvature fixed 
point, where all scales (except of course for the inverse 
ultraviolet cut-off) are larger compared to the radius of 
curvature. Performing the second derivative with respect 
to the conjugate field, we obtain the scaling expression 
for the order parameter susceptibility: 

x(t,K) = 

= <-'1-(^). (73) 

In the fiat space limit «: = 0, the scaling function behaves 
as <i>(a; —>■ 0 ) —)■ const, and we recover the fiat space 


results. On the other hand, our above analysis implies 
that for large argument $(a; ^ 1 ) cx holds with 
crossover exponent 


7 - 7 / 


24 


(74) 


Vf Ntt^ 

Here 7 is the susceptibility exponent of the hyperbolic 
space obtained above. Thus, we find %(/, k) cx . 

The behavior k~^ is, at the considered order, fully con¬ 
sistent with the (/)log(/c) behavior that occured in our ex¬ 
plicit analysis. We have calculated the critical exponents 
r], 7 and v at lowest order in 1 /N and found that these are 
identical to the exponents in fiat three-dimensional space 
at lowest order. For y and 7 we showed that 0(1/A^) cor¬ 
rections are absent. As our calculations show, the reason 
for this absence is the fact that correlations are exponen¬ 
tially decaying beyond the radius of curvature even at the 
critical point. The lowest order values of the exponents 
are computed from local quantities, which are oblivious 
to the finite curvature, whereas the higher-order correc¬ 
tions are determined through integration over the whole 
of hyperbolic space, wherein the finite curvature serves to 
cut off the long-wavelength fluctuations. For this reason, 
we may also surmise the absence of corrections to the 
other critical exponents. It is moreover plausible to as¬ 
sume for the same reason that higher-order corrections to 
the exponents will also be absent in the 1/iV-expansion. 
Thus we conjecture that the critical exponents we found 
are correct to all orders in 1/7V. 

An interesting question is how our results are modified 
for dimensions d different from 3. The Laplacian in d ^ 3 
dimensions is still gapped. The only modification in our 
lowest order 1 /A^ calculations of the critical exponents 
would be a change of the integration measure in (431, 
from to multiplied by a numerical factor. How¬ 

ever, this leads again the same saddle point equations as 
in fiat space. Thus we can make the stronger statement 
that all critical exponents in hyperbolic space are just the 
leading order 1/A^ exponents of fiat space. In particular, 
we have v = and 7 = for d < 4 and mean-field 
exponents for d > 4. The upper critical dimension is 
d = 4 even for finite k. 

In summary, we conclude that the description of many- 
particle systems in hyperbolic space is a promising avenue 
to investigate uniform frustration and non-trivial critical 
behavior within one theoretical approach. 
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